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I.  Introduction 


An  extrapolation  of  the  geomagnetic  secular  variation 
from  the  earth’s  surface  to  the  core-mantle  boundary  has 
previously  been  carried  out  for  an  insulating  mantle  (ref.l). 
An  improved  estimate  requires  a  knowledge  of  the  electrical 
conductivity  CT  throughout  the  mantle.  This  function  is 
roughly  known  only  to  about  the  700-km  depth.  With  greater 
depths,  it  is  presumed  that  CT  increases  monotonically.  The 
work  reported  on  here  is  an  investigation  of  the  radial 
distribution  of  O’  . 

The  work  has  been  divided  into  two  parts.  In  Part 
I,  we  employ  a  simplified  earth  model;  the  core  is  surrounded 
by  a  variably  conductive  medium  of  infinite  extent. 

Choosing  a  power  law  for  CT as  function  of  depth  the  equations 
are  explicitly  integrated  in  terms  of  Bessel  functions. 

A  surface  harmonic  expansion  at  the  core,  for  an  arbitrary 
frequency,  is  then  analytically  continued  to  a  distance 
equal  to  the  earth's  radius.  Neglecting  reflection  effects, 
we  thereby  obtain  the  wave  attenuation  and  phase  retardation 
after  propagation  through  the  mantle. 

Ourresults  show  that,  after  the  initial  rapid  rise 
in  the  first  700  or  800  km  of  the  mantle  (cf.  ref.  2),  the 
conductivity  levels  off,  G"at  the  1000  km  depth  is  of  the 
order  of  1/10  the  value  at  the  core-mantle  boundary. 

In  Part  II  we  consider  aperiodic  solutions  of  some 
relevant  models  together  with  the  continuations  from  Part  I. 
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II.  Mantle  Equations 

For  later  reference,  we  catalog  the  quasi-stationary 
field  equations  applicable  to  the  Earth’s  mantle.  The  general 
theory  of  electromagnetic  induction  in  non-uniform  conductors 
has  been  formulated  by  Lahiri  and  Price  (1939) .  The  following 
work  draws  to  some  extent  on  their  formalism. 

Equating  Maxwell’s  displacement  current  to  zero  we 
obtain  the  quasi-stationary  field  equations  (in  Gaussian  units) 
for  a  general  linear  isotropic  medium, 

V*  &=  -  \k  v  •  B  =  0  (1) 

Vx  H  =  ^  J  V  *  D  =  4ne  (2) 

The  material  equations  are 


D  =  KE,  B  =  UH,  J  =<tE  . 

From  the  first  of  Eqs.  (2),  it  follows  that  ^is  approximately 
solenoidal.  Neglecting  £  compared  to  S^K,  the  second  of  Eqs.  (2) 
gives  the  approximation 

4*6=  J  -V(-§>.  13) 

Thus  the  space-charge  density  vanishes  only  if  (a)  ^  lies 
in  the  iso-level  surfaces  K/<j=  constant,  or  (b)  if  K/(yis 
a  constant.  From  (1)  and  (2) and  mus t  satisfy  the  equations 


Vx  [  pVx  E, 
2 


+  kCjg,  =  0  (4a) 

.  +  kpH  =0  (4b) 

where  k  =  4u;/cc'.  On  the  other  hand,  in  terms  of  the  vector 
potential  we  have 

(5) 

where  A  is  any  solution  of 

Vx[^x  (pA)]  +  kpc?i  =  0  ,  (6a) 

subject  only  to  the  auxilary  condition 

V-  (pCT&)  =  o 

(Here  we  neglect  static  magnetic  fields.)  Our  choice  of  the 


(6b) 
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vector  potential  is  such  as  to  include  the  potential  of  the 
free  charge  distribution,  namely,  from  (2), 

-4ite=V-(f&)  •  (?) 

In  terms  of  we  list  below  the  boundary  conditions 
appropriate  to  the  transition  surface  separating  two  media 
M-^(K^,  Hi*  ^1^  and  ^2*  From  the  continuity  of  the  tangential 
components  of  JE,  and  E,  and  the  normal  component  of  J3,  we  have 


(8) 


$  x  *  *  *  {<PA>2  -  0 

*  *  A(p^pA)  =  0 

£  •  A  (\7x  pA,)  =0  . 

The  charge  dens ity  res iding  on  the  transition  surface 
may  be  expressed  in  terms  of  the  discontinuity  in  jD^,  as  well 
as  in  by  the  conservation  of  charge.  Thus  we  may  write 

=  [l— • 


where  n  is  directed  from  medium  1  into  medium  2. 
an  insulator  the  above  reduces  to 


(9) 

If  medium  2  is 


(10) 


Applying  these  equations  to  the  mantle  we  set  <T=  C Tip) , 
where^=  r/Rc  and  Rc  is  the  radius  of  the  core.  Hereafter  we 
shall  always  take|J  =  constant.  Under  these  conditions, 
standard  toroidal  solutions  have  the  same  form  as  in  the 
homogeneous  case  (0=  constant),  namely., 


A  =Vx 


x  r. 


*pcr4.  (11) 

General  poloidal  solutions  are  unknown.  These  solutions, 
however,  are  of  little  interest  in  the  study  of  the  secular 
variation  since  they  give  rise  to  toroidal  magnetic  fields 
which  are  not  observable  at  the  earth's  surface  (Elsasser,  1950), 
Standard  solutions  of  Eq.  (11)  are 

*.=  Rclfo</>’t>Vs"<e.0>  xr. 


(12) 


where  (^o,t)  satisfies  the  spherical  Bessel  differential 
equation,  with  variable  cr, 


a2 

bp 


i/flg)  +  kw** 


(13) 


S^(0,5/)is  the  usual  surface  harmonic. 

Outside  the  conductive  medium  (Jis  identically  zero, 
and  ife  is  computed  from  the  resulting  homogeneous  equation. 
Corresponding  solutions-^  are,  for  r^RQ  (earth* s  radius). 


i'nnT  Vn+1> 

•’f'nV 


(14) 


0  "1 

where  C  (t)  and  C  (t)  are  arbitrary  complex  functions  of  time, 

n#/ 

The  corresponding  vector  potentials  are,  for  r^R  , 


e 

^m 


n. 


k-~  “  R-I&t>frr)"l7s>,0)  x  r. 


n-1 


&L  =  VSn<e^  XX> 


(15) 


the  coefficients  being  chosen  so  as  to  satisfy  the  relations 
iL  =VX  .A,  =  Here  C^t)  represents  the  inducing  field 

whereas  C^t)  corresponds  to  the  induced  field  together  with 
the  field  of  internal  origin. 

For  later  reference,  the  components  of  A  and  H  are 


listed  below  (R  $r<R  ) 

C  G 


Ar  =  0 


R 

A  =  -  ? 

A0  slnB 


H 


,= /o  1  ^  n(n+l)\|4im 


n  ,m 


n  ,m 


Hq  =/<> 


(16) 


A*  =  -R. 


I 

n  ,m 


~  1<THh  Z  k>{P 


n,m 


Outside  the  mantle  it  is  more  convenient  to  replace  R  by  R 

C  6 

andyOby  r/Rg,  where  they  occur  explicitly  in  (16). 


III.  Periodic  Solutions  of  Mantle  Equations 


Aside  from  an  insulating  mantle,  the  simplest  earth 
model  is  that  of  an  infinite  conductive  medium  surrounding  the 
core.  We  wish  to  examine  this  model  with  the  chief  aim  of 
determining  roughly  the  distribution  of  the  bulk  of  electrical 
conductivity  in  the  mantle.  Waves  introduced  at  the  core 
boundary  are  investigated  as  to  their  attenuation  and  phase 
retardation  at  a  distance  equal  to  the  earth’s  radius,  reflection 
of  waves  being  neglected. 

For  the  electrical  conductivity  we  choose  the  three- 
parameter  distribution 

o{p)  =  c yf2  (17) 

Using  this  function,  we  make  the  following  substituions  in 
Eq.  (13),  (Compare  Lahiri  and  Price,  Eq.  5.3.) 

2  xi 


(18) 


=  i5/2(Wk^R c2c^Y 

U-  r/  * 

v„=  2[(n+i)2-^2q[/q)]j/tir-2i9  n  =  1,2,3,..., 

and  find  that  Wy  (z)  obeys  Bessel’s  differential  equation  of 

order  v„ and  argument  z.  For  simplicity,  in  these  initial 

calculations  we  shall  always  take  C£=0. 

In  the  computations  we  consider  only  the  radial  component 

H  for  a  sinusoidal  input  wave.  In  terms  of  the  constant 
r 

coefficients  if1  of  an  expansion  of  Hr  at  the  core,  we  then 
rn 

write,  from  (16), 


(19) 


where 


fW> 


"  (z) 


z  denotes  the  value  of  z  at  r=R  , 
c  c 


B  is  the  real  total 
n 


attenuation  and  0n  the  real  phase  retardation .  Since,  in  Part 
i,  a&>,  it  follows  that  v„  is  real.  Accordingly,  three  cases 
arise  depending  upon  whether  Tis  greater,  equal,  or  less  than  2. 


Case  I  (y>2).  The  equation  (18)  connecting^oand  z  shows  that^ 
->oo  implies  z-*-o .  Since  |A|=  0I/62  )  when  or  =0,  it  follows 
that  forC£>0,  must  vanish  to  at  least  this  order.  Therefore 
W(z)  must  vanish-  to  the  appropriate  order  when  z->o.  The 
only  solution  of  Bessel's  dif.  Bq.  satisfying  this  condition 
is  that  of  the  first  kind,  the  series  expansion  being  regular 
over  the  whole  cut  plane.  We,  therefore,  write,  for  large  values 


/=r^(/3't)  = 

#Hut0-3/8J 


(20) 


|^|y(mt+W«)  _n_2 

ru*»n)  |y-*r  “ 


=  0(/O 


r-n-2)  . 


Case  II  (ff=  2) .  This  case  is  singular  for  the  substitutions 
(18).  Eq.  (13),  however,  is  homogeneous,  and  a  direct  integration 
gives  the  general  solution 


=  °(1)  exp{1[^-bnlog/qJ- [an+3/2]lo^oJ+ 

( , 

C^2)  exp{i[o>t+bnlogyq|+  [an-3/s]  lo^oj  , 

where 

an  =^{s(n+3)2  +  ^{N4  +  (n+-s)4}}^n+i 

bn  =^/{-i(n+^)2  +  iV^4+(n+i)4]}^  °‘ 

We  are  interested  only  in  the  expanding  wavefront.  The  total 

attenuation  B  and  phase  retardation  0  in  Eq .  (19)  are 
n  n 

therefore 


K  -  =/53/2-a« 

For  large  J^j,  afi  i  bfi  =  |<£|2“*  so  that 

^n  ®  2~®|4?llogyO 

an  =(n+i>  *|£f(2n+l>-3, 


(22a) 


(22b) 


For  i«i  small. 


bn  *|£p/(2n+l)  and  therefore 


=l^|2(2n+1)"llo^> 

n  1  -n-2-  |£|4  (2n+l)" 
Bn  V3 


(22c) 


Case  III  (  2)  .  Referring  to  Eq.  (18),  we  see  that  with 

implies  |z|-*p  .  Since  the  attenuation  here  is 
greater  than  in  the  two  preceeding  cases,  it  follows  that 
Wv  must  vanish  as|z!->oo.  From  all  the  cylinder  functions  we 
must,  therefore ,  choose  one  of  the  two  Hankel  functions  since 
they  alone  vanish  in.  the  neighborhood  of  the  point  z  =  oo  . 
Consistent  with  our  choice  of g  and  the  positive  time  factor 
we  therefore  write,  in  terms  of  the  Hankel  function  of  the 
first  kind, 

=  ei“V'3/2H<li1)(z). 


For  large  |z| 


E{v1]  (z)/v[-|rei(2“w/2-7t/4) 


Therefore 


exp{-(l+i)2^|(/01-  y/2-i)/|r-2|}  . 


To  investigate  the  effects  of  a  material  medium  we 


define  the 


rsical  attenuation  C  as  the  ratio  of  the  total 
- n 


attenuation  13  to  the  geometrical  attenuation  G  : 

n  n 


n  -Bn 

’<5 


Gn  V° 


In  free  space  C  =  1*  Observe  that  in  the  asymptotic 
expressions  in  the  above  three  cases  Cn  is  (I)  constant,  (II) 
varies  as  an  invenoo  power  law,  (III)  varies  as  the  complex  of 
an  exponential  attenuation  and  a  diroet  power  law.  In  case  I, 
the  phase  retardation  is  constant;  the  conductivity  rapidly 
approaches  zero  and  ji  approaches  the  gradient  of  a  scalar  time 
dependent  potential,  the  equipotential  surfaces  varying 
isochronously. 

Tables  I  and  II  include  the  tabulated  values  of  Bn 
and  for  the  ^-values  indicated.  In  the  calculation  we  take 

5.4986*10B  [Syr/emu]  2  ,  U)  -  2tt/t 

6=C£/t  ,/S^E  C/Re  =  0.546875,  ^ 


TABLE  1.  TOTAL  AtTENUATION  Bn 

y  «  0,000 


6 

Bi 

B2 

B3 

B4 

B5 

B6 

B7 

B8 

B9 

B10 

B11 

B12 

0®6* 

,0103793 

,0103731 

— 

... 

— 

.0103216 

083 

.0B6425 

.O06291 

— 

— 

— 

.085188 

082 

.0*1634 

.061591 

... 

— 

— 

.0*1256 

081 

.0*1X05 

.041065 

... 

— 

.0*7607 

095 

.0*2168 

■  032058 

— 

... 

— 

.0*1273 

0®3 

•021091 

.021020 

... 

— 

... 

. 0*6467 

091 

.01128 

.01003 

... 

— 

— 

.023434 

0l06 

.02830 

.0241 

— 

— 

— 

.02574 

0101 

.08160 

,0664 

... 

... 

_ 

.02794 

-11 

0AA1 

.1623 

— - 

. — 

... 

— 

— 

0121 

.163 

.069 

— 

— 

— 

.  Q2  6  0 

Y  *  2.000 

061 

.0**2518 

... 

... 

— 

_ 

.0*2235T 

... 

_ 

... 

... 

071 

.0102621 

.0l02581 

•0102522 

,0102446 

,0102353 

.0102245 

.0102125 

.0101995 

.0101858 

.0101717 

.0101674 

,0101430 

O03 

.0*1051 

... 

— 

— 

... 

.0*7122 

_ 

_ _ 

_ _ 

.0*3377 

091 

.0*2367 

.032253 

•  032  09  0 

.031890 

,031664 

.0*1425 

.0*1186 

.0*9573 

.0*7501 

.0*5700 

.0*4203 

.0*3006 

093 

.026619 

,026039 

,025246 

.084325 

,023376 

. 022483 

.0*1744 

.021161 

.0*7405 

.0*4562 

.0*2734 

.0*1604 

091 

.03546 

•03007 

.02322 

.01631 

.01055 

. 026429 

.023765 

.022151 

.021211 

.0*6753 

.0*3745 

.0*2066 

0101 

.1421 

.08584 

.04615 

,02655 

.01457 

.027980 

.024368 

,022390 

.021307 

.0*7151 

.0*3911 

.0*2140 

O11! 

.1632 

.08941 

.04691 

.02675 

.01463 

.028000 

•024375 

. 022393 

.0*1308 

.0*7156 

.0*3913 

.0*2140 

0121 

.1635 

.08944 

.04891 

.02675 

.01463 

• 028000 

.024376 

,022393 

•021306 

.0*7156 

.0*3913 

.0*2140 

Y  *  9.440 


073 

,061186 

— 

0-- 

— 

— 

— 

.077067 

_ 

_ 

_ 

_ 

.072615 

072 

,052381 

— 

— 

— 

— 

— 

. 0*1161 

_ 

_ 

_ 

_ 

.0*3144 

07I 

,03U93 

,031069 

,049492 

047675 

. 046325 

.044923 

.0*3727 

.'0*2754 

.0*1991 

.0*1411 

.0*9818 

.0*6723 

0®7 

. 035488 

— 

... 

— 

— - 

.0*1332 

... 

_ 

--- 

_ 

.0*1960 

,083 

,027542 

,025962 

,024474 

,023220 

,022237 

.021510 

.0*9935 

.0*6396 

.0*4040 

.0*2510 

.0*1536 

.0*9281 

081 

.05138 

.03540 

.02322 

.01466 

,O20905 

.025376 

.0*3157 

.021827 

.021045 

.0*5928 

.0*3338 

.0*1869 

0®3 

.1298 

— 

--- 

— 

— 

“  Vact 

.024221 

... 

... 

_ 

... 

.0*2112 

0®1 

.1566 

.06766 

.04827 

.02650 

.01453 

2 

.0  3—0 

•024357 

.0*2385 

• 021305 

.0*7141 

.0*3906 

.0*2137 

0A01 

.1635 

.06942 

.04891 

.02675 

.01463 

,028000 

.024375 

.022393 

,021308 

.0*7156 

.0*3913 

.0*2140 

oni 

.1635 

.06944 

.04891 

.02675 

.01463 

. 028000 

.024375 

,022393 

•02130e 

.0*7156 

.0*3913 

.0*2140 

y  »  14.16 


,07S 

,07S 

.052212 

•044581 

•043604 

.042761 

,042069 

.0*1520 

.0*1097 

.0*4807 

■.037800 

.0*54  67 

.0*3784 

"■ 

.0*1752 

,07BB22 

.0B1174 

,071 

.022354 

.021653 

.0  2582 

.023270 

.021135 

.029802 

.0*507* 

.0*3316 

.0*2139 

.0*1363 

.0*8597 

.0*5368 

.0*3322 

:oa3 

.0*1 

.03740 

.02446 

.01565 

.026034 

.0*3660 

o 

,0*21B7 

.0*1297 

.0*7603 

.0*4417 

.0*2547 

.0*1463 

.1063 

.06453 

.03^6 

.0409? 

.02166 

.01231 

.0  6940 

.0*3887 

.0*2167 

.0*1203 

.0*6661 

.0*3680 

.0*2029 

,0*1 
,0101 
•  011! 

.1626 

.06906 

',02668 

-y  t 

.01460 

.0*7988 

.0*8000 

.0*4369 

.0*2390 

.0*1307 

.0*7150 

.0*3913 

.0*2139 

.1636 

.06944 

.04891 

.0264$ 

.01463 

.0*4375 

.0*2393 

.0*1308 

.0*7156 

.0*3913 

.0*2140 

.1635 

.06944 

.04691 

.02^5 

.01463 

.0*8000 

.0*4375 

.0*2393 

.0*1306 

.0*7156 

.0*3913 

.0*2140 

Y  -  18.90 

,0*1 

,051454 

,051079 

.067869 

*0*5686 

.064051 

. 0*2854 

.0*1990 

.0*1375 

.079421 

.  076404 

•074321 

. 072910 

,076 

.046845 

--- 

___ 

— 

--- 

— 

.0*7368 

— 

— 

— 

— 

■  0*8817 

,073 

. 036950 

.03486O 

.033279 

.0*2216 

,031465 

.0*9  634 

.0*6274  ' 

.0*4049 

.0*2593 

.0*1647 

.0*1040 

.0*6519 

,071 

.01384 

•O28904 

,025748 

im 

.022269 

.0*1403 

.0*8601 

.0*5229 

'  .0*3155 

.0*1891 

.0*1128 

.0*6666 

,ob5 

.04195 

.02616 

.01611 

.025890 

.0*3508 

.0*2072 

.0*1214 

.0*7065 

.0*4060 

.0*2349 

.0*1344 

.0*1 

.1393 

.07934 

.04400 

.02451 

.01359 

.0*7524 

.0*4153 

.0*2288 

.0*1259 

.0*6920 

.0*3797 

.0*2086 

,0*1 

.1632 

.06931 

.04666 

.02673 

.01462 

.0*7995 

.0*4373 

.0*2392 

.0*1308 

.0*7153 

.0*3912 

.0*2140 

,0101 

.1635 

.06944 

.04691 

.02675 

.01463 

.0*8000 

.0*4375 

.0*2393 

.0*1308 

.0*7156 

.0*3913 

.0*2140 

,0U1 

.1635 

.06944 

.04891 

.02675 

.014  63 

.0*8000 

.0*4376 

.0*2393 

.0*1308 

.0*7156 

.0*3913 

.0*2140 

y  •  20.33 


0®3 

.0*2588 

.0*17.49 

.0*1173 

.077809 

.075167 

,073398 

•072222 

.071446 

. 089360 

,08 6097 

.0*3860 

.082470 

0®2 

• 0344 12 

.0*2938 

.0*1942 

.0*1276 

.0*8328 

.0*5407 

.0*3491 

.0*2244 

.0*1436 

.079244 

,075802 

.073667 

0#1 

•031706 

,031106 

.047123 

. 044569 

.0*2900 

.0*1838 

.0*1156 

.0*7230 

.0*4611 

.0*2802 

.0*1733 

.0*1088 

07S 

.c^gooo 

,024966 

.0*3057 

.081871 

.0*1138 

.0*6891 

.0*4151 

.0*2489 

.0*1486 

.0*8838 

.0*5237 

.0*3099 

071 

.04936 

.02940 

.01739 

.01022 

.0*5977 

.0*3477 

.0*2014 

.0*1161 

.0*6670 

.0*3817 

.0*2177 

.0*1238 

0®3 

.1264 

.07142 

.04015 

.02248 

.01254 

.0*6878 

.0*3873 

.0*9146 

,0*1187 

.0*6383 

.0*3614 

..0*1991 

0®1 

•.1679 

.06661 

.04766 

.02616 

.01+3S 

.0*7867 

.0*4311 

.0*9381 

.0*1293 

.0*7062 

.0*3844 

.6*9192 

091 

.1636 

.06942 

.04890 

.02674 

.01463 

.0*7992 

.0*4374 

.0*2393 

.0*130. 

.0*7188 

.0*3(1! 

.0*2140 

0Ul 

.1636 

.06944 

.04891 

.02676 

.01463 

.0*8000 

.0*4376 

.0*9393 

.0*1308 

.0*7166 

.0*3913 

.0*2140 

•  Tbf  notation  «08S  danotaa  6  x  10"9 . 


TABLE  2.  PHASE  RETARDATION 

Y  *  O.OOQ 


0 

n6  * 

♦l 

22.773 

17.64 

14.395 

10.17 

7.180 

f2 

22.756 

17.61 

14.37 

♦3 

♦4 

*5 

*6 

22.61 

17.43 

14.15 

f7 

*0 

’9 

*10 

♦n 

*12 

u  0 

O03 

U  2 
O8! 

,095 

XO. 

7.132 

_ 

— 

— — 

6.73 

---- 

— 

.... 

.... 

.... 

i  O9  3 

5  491 

5.0 

mmmm 

_____ 

_  _  _  _ 

.... 

.... 

,091 

A10c 

5  550 

3.174 

2.216 

.9163 

.2013 

.03007 

74.16 

23.44 

3.06 

2.1 

.... 

,0  6 
n10i 

.0  1 
nll. 

.  o  1 
n12. 

iU  1 

A6! 

Y  =  2.000 

74.09 

.0  1 
.071 

23.43 

23.40 

23.37 

23.33 

23.29 

23.23 

23.17 

23.10 

23.03 

22.94 

22.85 

.081 

7.380 

7.34C 

7.267 

7.171 

7.054 

6.916 

6.760 

6.587 

6.400 

6.202 

5.996 

5.785 

.o9i 

2.259 

2.116 

1.924 

1.713 

1.508 

1.328 

1.176 

1.050 

.946 

.860 

.707 

.725 

o 

M 

O 

.525 

.355 

•  258 

.202 

.165 

.140 

.121 

.105 

.096 

.085 

.074 

.074 

.O11! 

.063 

.037 

.027 

.018 

.016 

.012 

.012 

.01 

.01 

.01 

.01 

.01 

.O12! 

.006 

.00 

.00 

.00 

.00 

.00 

.00 

.00 

.00 

.00 

.00 

.00 

Y  =  9.440 


,073 

16.22 

.... 

— 

.... 

— 

15.58 

— 

— 

.... 

— 

14.81 

,071 

9.310 

9.176 

9 . 126 

8.863 

8.700 

8.535 

8.369 

8.203 

8.030 

7.875 

7.715 

7.556 

0e7 

7.729 

— — 

— 

— 

6.742 

— 

.... 

— 

— 

5.947 

,083 

4.880 

4.704 

4.529 

4.357 

4.189 

4.026 

3.869 

3.718 

3.572 

3.433 

3.300 

3.172 

,031 

2.545 

2.366 

2.197 

2.040 

1.896 

1.765 

1.646 

1.538 

1.441 

1.354 

1.275 

1.204 

,093 

1.055 

— 

— 

— 

— 

.5321 

— 

— - - 

---- 

— 

.3706 

.091 

.3802 

.3221 

.2770 

.2436 

.2175 

.19  63 

.1789 

.1643 

.1519 

.1412 

.1320 

.1238 

,o10i 

.0385 

.0323 

..0278 

.0244 

.0218 

.0197 

.0179 

.0164 

.'0153 

.0141 

.0132 

.0124 

.d11! 

.02389 

.02327 

•  022'81 

.02247 

.02220 

.02199 

.  02181 

.02166 

.02154  * 

. 02143 

.02134 

,02125 

,012! 

.03389 

.03327 

.0S281 

— 

.... 

— 

— 

— 

.... 

—  - 

— 

•03125 

Y  *  14.16 


*$» 

0*3 

rt 

13.66 

..... 

— 

... 

... 

— 

12.93 

_  ___ 

_ 1 

____ 

12.34 

10.45 

10.33 

10.21 

10.09 

9.966 

9.846 

9.727 

9,609 

9.493 

9.377 

9.263 

9.150 

072 

8.430 

— 

— 

— 

M 

— 

_ 

_ 

____ 

____ 

_ _ 

071 

083 

5.790 

5.666 

5.545 

5.426 

5.194 

5.081 

4.971 

4.862 

4.756 

4.653 

4.551 

2.098 

2.778 

2.663 

2.551 

2.444 

2.341 

2.245 

2.150 

2.061 

1.977 

1.897 

1.821 

081 

1.380 

1.277 

1.178 

1.090 

1.012 

0.9452 

.8818 

.8271 

.7020 

.7344 

.6950 

.6594 

091 

0101 

.1635 

.1446 

.1295 

.1173 

.1072 

.09871 

.09145 

.08518 

.07926 

.07492 

. 07060 

.06681 

.01639 

.01448 

. 01297 
.021297 

.01173 

.01074 

. 029882 

.029155 

.028527 

. 027979 

.027498 

.027072 

.026691 

•  011! 

.021640 

.021448 

.021174 

.021074 

.039882 

Y  =  18.90 

•039155 

. 038527 

. 037979 

. 037498 

. 037072 

.036691 

.061 

14.00 

13.92 

13.83 

13.74  . 

13.65 

13.56 

13.47 

13.38 

13.29 

13.20 

13.12 

13.03/ 

.071 

4.055 

3.966 

3.878 

3.795 

2.46? 

3.706 

3.625 

3.544 

3.465 

3.387 

3.311 

3.237 

3.164 

•  085 

2.703 

2.615 

2.529 

2.365 

2.286 

2.211 

2.138 

2.067 

2.013 

1.934 

1.072 

.081 

.8407 

.7660 

.7104 

.6684* 

.6232 

.5855 

.5509 

.5200 

.4922 

.4  671* 

.4444 

.4237 

.  091 

.0898 

.0816 

.0748 

.0690 

.0641 

.0598 

.0561 

.0528 

.0498 

.0472 

.0448 

.0426 

.0101 

. 02899 

.02817 

.02749 

.02690 

.02640 

,02598 

. 02561 

.02528 

.02499 

.02472 

.0Z448 

.02427 

.011! 

,0390 

,03S2 

.0375 

«0369 

.0364 

.0360 

.0356 

.  0353 

♦  0350 

.0347 

.0345 

.0343 

.0121 

.0*9 

.0*8 

.0*7 

.0*7 

.0*6 

.0*6 

Y  =  28.33 

.0*6 

.0*5 

.0*5 

.0*5 

.0*4 

.0*4 

.  063 

15.70 

15.64 

15.58 

15.52 

15.46 

15.40 

15.34 

15.28 

15.21 

15.15 

15.09 

15.03 

.062 

12.72 

12.66 

12.60 

12.55 

12  .,49 

12.43 

12.37 

12.31 

12.25 

12.20 

12.14 

12.00 

.061 

8.851 

8.792 

8.732 

8.675 

8.617 

8.559 

8.502 

8.445 

8.388 

8.332 

8.276 

8.219 

,073 

4.621 

4.562 

4.505 

4.447 

4.391 

4.335 

4.280 

4.225 

4.171 

4.117 

4.065 

4.013 

.071 

2.449 

2.359 

2.335 

2.279 

2.225 

2.172 

2.119 

2.069 

2.019 

1.971 

1.-923 

1.870 

.083 

1.053 

1.002 

.9542 

.9102 

.8695 

.8310 

.7968 

.7643 

.7341 

.7059 

.6797 

.6552 

,081 

.3857 

.3620 

.3410 

.3222 

.3054 

.2901 

.2722 

.2638 

.2523 

.2417 

.2301 

.2231 

.091 

.03914 

.03665 

.03445 

.03250 

.03075 

.02919 

.02778 

.02650 

.02533 

.02426 

.02328 

.02237 

.011! 

•0S3»16 

. 033666 

.033446 

•033251 

,033076 

. 032920 

•032779 

. 032  650 

.032534 

.032426 

•032328 

.032237 

*  The  notstlon  .p86  denotes  5  x  10  B. 


where  <r  is  measured  in  emu  and  T  in  years.  is  dimensionless. 

^  (1) 

Standard  asymptotic  expansions  of  H  and  were  employed  for 

those  larger  values  of  6  for  which  the  series  e xpansions  failed 

to  converge  rapidly  (approximately  15  terms).  The  values  of 

Bn  and  are,  in  general,  accurate  to  four  significant  figures. 

Computational  errors  were  checked  only  by  continuity  of  graphs. 

The  values  of 2T=  9.440,  14.16,  18.90,  28.33'were  chosen  such 

-2  -  3 

that  the  conductivity  at  the  700-km  depth  would  be  10  ,10  , 

-4  -6 

10  ,  10  times  O' ,  the  value  at  the  core  surface. 

’  o’ 

Figures  1  to  5  show  plots  of  the  first  twelve  'coefficients.. 

over  the  appropriate  range  6,  or  ,y shown.  The  left-hand 

horizontal  portion  of  each  curve  indicates  the  geometrical 

attenuation.  Observe  that  with  increasing  'S’,  the  curves  are 

shifted  to  the  right.  That  is,  for  a  constant  period  of  vibration 

and  total  attenuation,  an  increase  in  ^"requires  an  increase  in 

CTq.  The  fact  that  the  first  and  twelfth  harmonics  approach 

each  other,  as  6  increases ,  clearly  shows  that  C  increases 

n 

with  increasing  order  of  harmonic - -lower  order  harmonics  are  physic*.lly 

attenuated  more  than  those  of  higher  order. 

In  Fig.  6  the  total  attenuation  of  the  first  harmonic 
is  plotted  as  a  function  of  (Tthroughout  the  range  0  to  35. 

Fig.  7  shows  the  physical  attenuation  for  the  first, 
seventh,  and  twelfth  harmonics  for  values  of  2.000,  9.440, 
and  28.33.  With  lncreas ing  Breach  family  becomes  more c  ompact 
and  is  displaced  towards  the  right. 

Fig.  8  shows  the  relative  physical  attenuations 
Cy/C-^  and  C-^g/C^.  For  example,  for  9.440  and  6=  10“8emu/yr., 
we  see  that  the  first  harmonic  is  attenuated  approximately 
44  times  the  twelfth  harmonic. 

The  phase  retardation  0n  is  plotted  in  Fig.  12  for 
^-values  equal  to  2.000,  9.440,  18.90  and  28.33.  The  travel 
time  required  for  waves  leaving  the  core  to  reach  the  earth's 
surface  is  obtained  from  the  expression 

t  =  0nT/2Tt,  (26) 

where  T  is  the  wave  period.  For  example,  for  2f=  9.440  and 

Q 

8  =  10  emu/yr.,  we  find  for  the  first  and  twelfth  harmonics, 
t^  =  1.48  yrs.,  t^g  =  1.20  yrs.  so  that  the  first 
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harmonic  would  arrive  approximately  3  l/3  months  after  the 
twelfth  harmonic.  This  quasi- dispersive  effect  occurs  for  a  pure 
sinusoidal  source. 

IV,  Equivalent  Conduc tlvlty 

The  utility  of  an  equivalent  § (equivalent  conductivity) 
lies  in  the  fact  that  calculations  involving  the  earth's  mantle 
may  be  simplified  by  a  choice  of  the  radial  dependence  of  the 
electrical  conductivity,  such  as  a  constant  or  an  inverse  square 
law.  We  next  establish  such  an  equivalence  for  later  reference. 

For  a  harmonic  of  given  degree.  Fig.  9  shows  that  curves 
Bn,  for  different  values  of  £f,  are  nearly  "parallel**.  This 
implies  that  one  may  establish,  for  each  harmonic .  an  equivalent 
S»  6c,  valid  over  a  reasonably  large  Bn~interval  by  a  horizontal 
translation  of  that  Bn-curve  for  which  £==0  (constant  electrical 
conductivity).  The  translations  log^'T  =  log1Qg  -  log105c  of  the 
zero-curve  to  any  chosen  ZT-curve,  for  various  constant  physical 
attenuations,  are  shown  in  Fig.  10.  Curves  I,  II,  III,  for 
n  =  6,  correspond  to  constant  physical  attenuations  0.800, 

0.500,  and  0.204.  The  attenuations  of  the  first  harmonic, 
corresponding  to  these  constant  attenuations,  are  shown  in  Fig. 

9  by  the  curves  I' ,  II' ,  III' ,  and  the  geometric  means  between 
the  extreme  points  (  0  and  ZT  =  35)  of  these  curves  are 

shown  by  curves  I",  II",  III".  Here  the  physical  attenuations 
are  (J-^  =  0.37,  0.18  and  0.061.  The  displacements  log  6/6c 
are  shown  In  Fig.  10,  for  n  =  1,  for  the  curves  I**,  II**,  III**, 

together  with  the  curve  III.  Curve  III  represents  a  total 

-3 

attenuation  of  amount  Bn  =  1.63  x  10“  ,  the  corresponding 
physical  attenuations  being  =  0.010  and  Cg  =  0.204.  In 
Fig.  10,  observe  that  the  curves  for  n  =  6  are  g rouped  closer 
together  than  those  for  n  =  1.  Curves  for  n  greater  (lesser) 
than  6  lie  below  (above)  those  for  n  =  6.  An  Inversion  of  the 
curves,  n  =  6,  occurs  at  8.5,  curves  I  and  III  intersecting 
Curve  II.  The  results  for  n  =  12  have  not  been  inducted  here. 
However,  from  the  graphs  of  the  physical  attenuation.  Fig.  7, 
one  may  expect  this  family  to  be  quite  compact  and  to  lie 
somewhat  closer  to  the  family  n  =  6  than  does  the  family  n  =  1. 
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Data  scaled  from  Fig.  9  is  given  below.  The  tabulated  values 
indicate  l0gj_Qrf  •  The  last  three  columns  pertain  to  the  family 


I 

II 

TABLE  3 

III 

I" 

II" 

III" 

2.000 

0.407 

0.38 

0.337 

0.435 

0.38 

0.35 

5.000 

---- 

1.04 

0.93 

0.84 

9.440 

1.25 

1.24 

1.22 

1.61 

1.48 

1.38 

14.16 

1.50 

1.50 

1.51 

1.98 

1.85 

1.76 

18.90 

1.68 

1.70 

1.71 

2.25 

2.12 

2.03 

28.33 

1.99 

2.00 

2.01 

2.59 

2.47 

2.14 

35.00 

_  _  _  _ 

—  a* 

^ 

2.78 

2.68 

2.59 

In  general,  different  attenuation  curves  cannot  be 
made  to  coincide  everywhere  by  a  relative  “horizontal  translation. 
An  attenuation  curve  Bn(§;^)  coincides  with  the  translated  zero 
curve  Bn(T  6^0)  only  in  the  ^-neighborhood  about  the  point  of 
intersection  of  these  curves  with  the  specified  physical 
attenuation  which  characterizes  T'.  To  indicate  the  error 
incurred  by  using  the  translated  zero-curve.  Fig.  lib  shows, 
roughly,  the  maximum  and  minimum  ratio  of  Bn  ( Sj^T)  to  Bj 
taken  along  the  same  8-ordinate  for  both  n  =  1  and  n  =  6, 
over  the  range  of  attenuation  in  Fig.  9.  Throughout  the 
Tf-range  indicated,  the  maximum  error  occurs  at  approximately 
the  same  value  of  total  attenuation.  These  values  are 
indicated  by  the  parenthesized  quantities  in  Fig.  lib. 
corresponds  to  curve  III. 

Fig.  11a  again  illustrates  the  maximum  and  minimum 
values  of  the  ratio  Bn  ( S * )  to  (<f §c ;0 )  for  the  constant 
attenuation  shown  by  curve  III.  Here,  however,  the  roles  for 
n  =  1  and  n  =  6  were  interchanged  in  computing  ff'.  Thus  the 
curves  marked  n  =  1  (n  =  6)  in  Fig.  11a  indicate  the  maximum 
error  in  the  first  (sixth)  harmonic  due  to  the  employ  of  the 
tranp lated  8  computed  from  the  sixth  (first)  harmonic. 

We  postulate  that  the  secular  field,  Hr(Re,6,0,t) 
would  require  twelve  harmonica  for  a  representation  with  an 
error  not  larger  than,  say,  3  or  4$.  The  translation  of  an  inter- 
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mediary  harmonic,  say  n  =  6,  would  therefore  be  expected  to 

represent,  very  roughly,  the  translation  of  the  representation 

of  the  whole  secular  field.  Since  the  'T'-curves  for  n  =  6, 

Pig.  10,  are  quite  compact,  we  arbitrarily  choose  curve  II  in 

the  approximation  of  the  equivalent  S,  &>•  From  the  Maclaurin 

2  3  ^ 

expansion  log  ,  nY  =  a  o  +b7  +cu  +. .. ,  we  may  write 
-LU  T 

C T=  Tc5c  =  exp  £-(aV+blf2+cft3+...)logeloJ,  (27) 

where  a,b,c,...,  are  determined  from  Table  3.  For  curves  II 
and  IIW  the  first  few  coefficients  are  tabulated  below.  For 
these  curves,  the  points  2.000,  9.440,  18.90,  28.33  and 
2.000,  14.16,  28.33  were  employed. 


CURVE  II 
a  0.211 
b  -0.0113 
c  0.000324 


CURVE  II** 
0.202 
-=0  •  00399 
6.89  x  10”5 


d  -3.60  x  10~6  - 

Of  especial  interest  is  the  value  Tf=  5.45(cf.  Eq.38).  From 

(27),  with  T  =  T,  we  compute  the  equivalent  conductivity  in 
c 

terms  of  the  conductivity  C£  at  the  core-mantle  boundary 

CT  =  0.1370  (JQ ,  Curve  II(n  =  6)  ( 

=  0.1164  ,  Curve  II”  (n  =  1) 


In  Part  II  we  shall  make  use  of  (28)  to  compute  from  a 
direct  computation  of  C£. 

V.  Computation  of  Electrical  Conductivity 


If  one  neglects  reflections  at  the  earth's  surface, 
the  electrical  conductivity  of  the  mantle  may  be  estimated  by 
a  method  which  replaces  the  time  averaged  squared  field  of 
at  the  core  by  a  properly  weighted  static  distribution.  The 
method  has  the  advantage  that  a  knowledge  of  the  absolute 
magnitude  of  the  field  at  the  core  is  unnecessary;  only  the 
relative  field  is  required.  We  first  consider  a  general 
expansion  of  Hr  in  the  form  ^  (29) 

H r(/>,  6,0, t)  =  ^~  Bne”i^n  ^  P*(cose)Lnmcosm0+bnmsinm0 
n^T  raSCT  L 


9 
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where,  from  Eq.(lQ,  a  and  b  are  proportional  tol/Cl.t)  and 
'  ^  nm  nm  i* 


Bn  and  0n  are  again  the  real  total  attenuation  and  real  phase 
retardation.  Making  use  of  the  usual  orthogonality  relations, 
the  surface  average  of  JhJ  is 


n=l  L 


£t  . 


Observational  values  of  Bn  (Re/Rc  ,t)  anm  and  Bnbnm  are  provided 


by  Vestine  (ref.  6, 


In  terms  of  associated  Legendre 


polynomials,  his  analysis  is  written 


Pnm<-e)[A”cos^+  (31) 

n=l  L  m=l  L  J  J 

where  Am  and  Bm  are  tabulated.  According  to  (31)  the  time 
n  2  n 

average  of  h  is 


r/LZn 


nt 

Zn+I 


V  =  (2li) 

n-fc  n 


The  subscript  t  denotes  the  time  average. is  tabulated 

below.  Table  4,  for  n  =  1,2 ,...,6.  The  time  average  was  taken 

only  over  the  four  ten  year  epochs  1912.5,  1922.5,  1932.5  and 

1942.5.  The  values  of  Am  or  Bm  found  separately  from 

n  n 


vj  tf/yi 

46.58 

89.67 

69.44 


TABLE  4 


Vn|?/yi 

51.01 

56.17 

33.31 


Vestine »s  X  and  Y  analyses  were  averaged  algebraically.  The 

significance  of  some  of  these  coefficients  is  somewhat  in  doubt 

because  of  the  rather  large  discrepencies  found  in  his  two  values 

for  the  same  coefficient.  The  error  in  V  f  is  further  augmented 

n  t 

by  the  limited  data$  the  epochs  for  which  there  exist  reliable 
date  cover  only  40  years.  To  supplement  Table  4,  a  straight 
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line  (y  =  -12.26n  +  110.3)  was  fitted  by  the  method  of  least 

squares  to  the  points  (Vn^,n)  for  n  =  2, 3, 4,5, 6.  These  results, 

together  with  the  first  term  of  Table  4,  are  shown  in  Table  5. 

A  A 

Plots  of  ,  for  each  epoch,  and  Vn^  are  shown  in  Pig.  15a. 


TABLE  5 


n 

1 

2 

3 


y(2?yr.) 

46.58 

85.80 

73.54 


n 

4 

5 

6 


y  (7f/yr . ) 

61.28 

49.02 

36.76 


Evidently,  for  the  particular  time  interval  represented  by 
these  plots,  the  activity  of  the  secular  field  successively 
diminishes . 

O 

We  next  compare  (32)  with  the  time  average  of  h  in  (30). 
Since  the  two  expressions  are  identical,  term  by  term,  we  then 
have,  at  the  earth's  surface. 


VVnt 


nm 


It* 


lu  ,  2*1  (n+m)0‘z 
PnmltjTHTSnj 


(33) 


At  the  core,  Bn  =  1,  ^  =  0  so 


that  Hr  is  expressed  in  terms  of 


a_  and  b  •  Various  criteria  are  available  for  estimating 
nm  nm  ° 


the  quantities  a 


nm 


These  will  be  discussed  shortly.  The 

,th  .th 


nnr 

ratio  b-j/bj  °f  the  total  attenuation  of  the  iUii  to  the  solid 

harmonic  is  now  compared  with  the  same  ratio  calculated  from 

Table  I.  Figures  13,  14  and  15  show  all  these  calculated 

ratios  for  n  =  1,2,. ..,6.  Also  included  is  the  ratio  B,0/B,. 

. ,  w  Vauue  of  .  id  l 

In  these  graphs,  a  specified  ratie  B^/B^  intersects  the  various 

2f-curves  at  critical  values  of  §.  Curves  for  these  points  on 

a  8  coordinate  system  for  various  number  p airs  (i,j)  may  be 

expected  to  intersect  in  a  small  region  of  the  plane  (ideally, 

at  a  single  point  if  only  one  frequency  were  present;  increasing 

the  frequency  would  Increase  the  ordinate  of  this  point.)  The 

values  of  ?T and  8 appropriate  to  the  mantle  are  then  read  off  the 

graph. 

We  now  return  to  the  coefficients  a  and  b  .  Because 

nm  nm 

of  the  slow  convergence  of  (29)  at  the  core,  we  choose  as  a  first^ 
rough  approximation 

u  ,2|2  (i+j)  i  I  _  ,2  |  2  (k+Jt)i| 

l*i j  It  jsm  TTrjyrJ-  halt  pri  TE^cnJ  * 
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for  all  i,j,k,X.  Accordingly,  the  ratios  B^/B.  In  Eq.  (33)  are 


L<1+1>  v3tJ 


Hz 


(34) 


Further  refinements  in  these  ratios  will  now  be  obtained  from 
a  consideration  of  the  magnetic  field  at  the  core. 

The  bulk  of  the  magnetic  topography,  for  Hr,  has  been 
described  in  terms  of  an  extensive  system  of  ridges,  in  motion, 
the  function  Hr  varying  rapidly  as  we  traverse  the  core  surface 
(ref.l).  The  cross-sectional  shape  of  a  typ-leal  ridge  at 
present  is  unknown.  The  ridge  half -width,  however,  is  sufficiently 
small  (3  to  10  degrees  of  arc)  to  make  the  expansion  (29) , 
at  the  core,  quite  insensitive  to  the  cross-sectional  shape. 

To  illustrate  this  insensitivity  we  tabulate  the  first  few 
coefficients  of  a  zonal  harmonic  expansion  of  three  simple 
functions  of  axial  symmetry. 


f,  = 


f„  = 


A, 


9 


A 


10 

The 


CO30-COS°< 
1  -  COSo< 

-  »  » 

fx  -  0,  c<  $  e<  % 

l ,  e^c* 

i  f2  =  o,  o<<e<7i 

1-CO30 
1-C03O4  * 

6<o<  5  f5  = 

o  ,  o< <e< % 

f  i«=15° 

o 

o 

CO 

il 

fpfc=15°) 

f  p(<*=30°) 

f,  (<*=30° ) 

X 

xO  .500 

x0.500 

o 

0.00852 

0.0335 

0.00852 

0.0335 

0.0335 

0.0253 

0.0960 

0.0251 

0.0939 

0.0915 

0.0412 

0.146 

0.0404 

0.135 

0.125 

0.0556 

0.177 

0.05  37 

0.150 

0.124 

0.0683 

0.186 

0.0641 

0.137 

0.0881 

0.0787 

0.174 

0.0715 

0.0993 

0.0251 

0.0865 

0.144 

0.0755 

0.0467 

-0.0502 

0.0916 

0.102 

0.0755 

-0.0088 

-0.120 

0.0938 

0.0564 

~  —  w.  — 

0.0932 

0.0138 

0.0900 

-0.0197 

examples  serve  to  illustrate  the  increase  of  insensitivity 


to  cross-sectional  shape  with  decreasing  half-angle.  For 
functions  of  half --angle  5°-10°  the  effect  is  considerably  enhanced. 
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Neglecting  the  cross-sectional  shapes  of  the  ridges,  it 

I  ^  i  /  ry  n  \  -  _  _  .  T _ *  ' 


la  I2  and 

|b  | 

1  nmlt 

1  nm| 

approximated  by  thef corresponding^  (^efficients  computed  from  a 
properly  weighted  static  distribution  which  varies  rapidly  over  the 
core.  For  a  purely  random  distribution  of  the  function 


Hr(Rc,S)  =  An?n(cos9)  (35) 

n 

of  axial  symmetry,  we  make  use  of  the  Legendre  addition 
theorem  (Eq.  4,  ref.  1)  and  compute  the  ratio 

Bj/Bj  =  lA/A^  •(Vit/Vj^)*  (36) 


which,  again,  is  independent  of  the  magnitude  of  the  field  at 
the  core.  For  small  half-angles  |!U|  =  |Aj|  so  that  (36)  does 
not  differ  appreciably  from  (34). 

Our  choice  of  a  random  distribution  seems  quite  justified. 
Vestine's  secular  variation  maps  show  a  relatively  rapid  regional 
motion  of  the  ridges  superimposed  upon  the  general  westerly 
drift.  Taken  over  a  few  hundred  years  one  would  therefore  not 
expect  a  0-  dependence  in  the  time  average  of  |h  |  2.  Like¬ 
wise,  there  is,  at  most,  a  small  ©-dependence ;  the  4  or  5 
active  localized  regions  at  the  core  occur,  evidently,  quite 
randomly.  £If,  however,  these  regions  remain  within  the 
latitudes  ,  say  75°N  and  75°S,  the  field  at  the  earth’s 
surface  would  still  be  large  over  the  polar  caps  and  would 
therefore  tend  to  mask  any  small  ©-dependence.]  The  ©-dependence 
at  the  earth’s  surface  may  be  illustrated  by  the  line  integral 
of  Mt  along  the  parallels  of  latitudes 


Table  6  below  shows  f(©)  computed  from  Vestine’s  data 
(ref.  5,  pp.  304-311)  for  10°  intervals  of  latitude.  The 
time  average  includes  only  the  four  available  epochs,  which 
are  insufficient  to  depict  the  expected  equatorial  symmetry. 

For  a  comparison  with  the  random  distribution,  Eq.(36), 
the  following  interesting  model  illustrates  a  O-dependence: 

Two  identical  functions  (35)  of  axial  symmetry  and  half -angle  oc 


-> 


are  placed  adjacent  to  each  other  at  the  equator  in  a  fixed 
east-west  orientation,  their  centers  being  separated  by  an 
arc  2  <*.  one  function  is  then  reflected  at  the  core  surface 
(An->-An)  such  that  their  common  point  is  a  point  of  reflection. 
i  Denote  the  (6,00  coordinates  of  the  centers  of  these  functions 

by  the  subscripts  1  and  2.  This  "dipole*  is  then  distributed 
over  the  core  surface  according  to  the  restraining  equations 

e1  =  e2  =  e,  =  0  +  of,  02  =  0  -  . 

At  the  poles  the  centers  coincide  and  the  functions  just  cancel. 


< 

e 

f(6)  ^yr. 

TABLE  6 

e 

f(6)  V/jr. 

0 

39.5 

180 

53.4 

10 

39.7 

170 

68.9 

\ 

20 

40.8 

160 

88.8 

30 

46.3 

150 

92.9 

40 

56.2 

140 

81.2 

V 

50 

59.1 

130 

65.3 

60 

55.8 

120 

54.6 

70 

53.1 

110 

55.2 

80 

56.8 

100 

60.0 

90 

60.9 

1 

■  3  ‘ 

Using 

total 

the  Legendre  addition 

attenuation 

theorem  we 

compute 

the  ratios  of 

if 

VBj =  hj  • 

lVAi|-  <v 

n/V* 

(37) 

9 

where 

l 

jr 

*  _r(2i+D. 

i  -  prsD 

(2j+l)sino< 

m+lTsTmx 

-sin (2 j+] 
-s in  121+] 

£]* 

If  were  unity,  (37)  would  be  equal  to  (36). 

Pigs.  15b  to  15i  show  the  plots  of  the  ratios  B^/Bj 
computed  from  Eqs.('34),  (36),  and  (3rI),  using  Tables  4  and  5.  The 
quantities  An  in  (36)  and  (37)  were  taken  from  Appendix  A, 
ref.  1,  for  the  appropriate  half-angles  vindicated  in  the 
figures.  The  figures  pertaining  to  (37)  have  been  included 
only  for  their  heuristic  valuej  the  static  distribution  at 
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the  core  is  too  unrealistic .  In  all  the  figures,  except  15d  and 
15e,  the  region  of  greatest  convergence  of  the  curves  was  roughly 
estimated  by  locating  graphically  the  minimum  of  the  relative  mean 
deviations  of  the  points  of  intersection  of  the  curves  with  the 
3~-ordinates .  The  relative  mean  deviations  were  computed  for 
increments  of  7T =1.25.  The  average  critical-values  of  5*and  HT 
are  (1  ohm”^  m”  ^=10’"^' emu) 

=  5.45  R  =  7.20  x  10_10emu/yr.  (38) 

Here  the  data  of  Pig.  15i  was  weighted  by  a  factor  of  2.  If  it 
is  assumed  that  Vestine’s  coefficients  Am  and  Bm  are  representative 
of  an  average  period  of  50  years,  the  electrical  conductivity  at 
the  core  -mantle  boundary  and  the  equivalent  conductivity  from 
Eq.  (28),  are  (39) 

CT  =50  yr  x  £  =3.60  x  10~8emu,  CJ  =4.93  x  10“9emu. 

o  ^av o  *  c 

-9 

A  lower  limit  x  -10  emu;-  obtain od  by  -oho»e-iHf»g'-an 

average  peit*a.»d  equa4— to  -the . epoah-  ir.tfiti.val,  namely..,  ■lfl~.yr<g-. 

According  to  the  assumed  distribution  of  conductivity 
the  value  of  ?f  in  Eq.  (38)  implies  that  the  conductivity 

at  the  core-mantle  boundary  is  only  10.6  times  that  at  the  -1,000  km 

-9 

depth  10  omu)  ;  below  this  depth  (Tis  essentially  uniform. 

At  lesser  depths,  one  would  hardly  expect  the  secular  field  to 
be  successfully  employed  in  the  computation  of  CTsince  the  physical 
attenuations  are  scarcely  influenced  by  the  outer  1/3  portion 
of  the  mantle.  At  the  earthfs  surface  typical  values  of  CTare 
as  follows:  wet  ground,  10~^  to  10"^^emuj  dry  ground,  10"*^  to 
10  ^erau;  fresh  water,  10“^emu$  sea  water,  4  x  10”^^emu. 

Compatible  with  the  results  of  Lahiri  and  Price  (1939),  it 
follows  thatCmust  increase  very  steeply  in  the  outer  1/3  portion 
of  the  mantle. 

The  tentative  results  (38)  and  (39)  are  considered  here 
as  only  orders  of  magnitude;  the  precision  has  yet  to  be 
established.  This  will  entail  an  Investigation  of  the  various 
factors  influencing  the  critical  values  of  Tf  and  6  •  Unfortunately, 
several  of  the  computed  ratios  B^/B y  for  which  i-j,  are  in¬ 
compatible  with  the  ratios  computed  from  Chapter  III.  These  ratios 
are  indicated  by  an  asterisk  in  Appendix  D.  Some  difficulty 
lies  in  the  scanty  data.  Of  particular  trouble  is  the  small 

dipole  term  which  may  be  indicative  of  a  weak  9-dependence( Observe 
that  no  difficulty  exists  for  the  ratio  Bg/B-^  when  (37)  is 
employed. ) 
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The  final  conclusions  of  this  chapter  will  be  contained 
in  Part  II  of  this  report. 

VI.  Observatory  Records 

Further  information  concerning  the  nu.-'re  of  the 
electromagnetic  waves  passing  through  the  mantle  may  be  found  from 
an  examination  of  the  magnetic  observatory  records.  Vestine's 
maps  of  the  secular  variation  fields • represent  10  yr.  averages 
centered  at  the  various  epochs.  Since  many  observatory  records 
show  disturbances  of  much  shorter  duration,  considerable 
activity  is  not  represented  in  these  maps.  Figures  16  through 
22  illustrate  this  point.  The  values  of  H  (horizontal 
component)  and  H  (=  -2)  were  scaled  from  Ves tine * s  graphs  of  the 
time  variation  of  the  earthds  main  field  (ref.  4,  pp.  151-249). 

The  solid-  and  open-circles  in  Figs.  16,  18,  20,  22  indicate 
map  values  of  fi  and  ,  respectively  .  Evidently  a  rapidly 
moving  ridge  at  the  core  would  not  be  recorded  on  the  maps. 

The  examples  clearly  show  that  no  appreciable  relative 
amplitude  attenuation  occurs  for  wave  periocb  greater  than  three 
or  four  months . 

All  the  stations  represented  by  Vestine’s  data,  99 
magnetic  observatories,  were  examined  for  a  possible  phase 
relation  between  and  ft  components.  Of  the  stations 
represented,  16  were  deleted  because  of  insufficient  data 
(either  the  H-  or  Z-  component  was  missing  or  the  time  interval 
covered  was  too  short.)  The  remaining  stations  were 
catalogued  according  to  the  phase  relation  between  the  vertical 
and  horizontal  components:  25  stations  were  characterized 
strongly  by  a  ”180°  phase  difference";  Hr  and  -H  increase 
and  decrease  together.  Figures  16,  17  and  22  are  typical 
examples.  Twenty  additional  stations  were  counted  as  being 
typified  by  the  180°  phase  difference  over  a  large  portion  of 
the  range  1905  to  1945,  but  with  exceptions.  Figures  18  and 
19  are  examples  of  this  latter  class,  although  several  of  the 
observatory  records  counted  here  were  not  as  consonant. 

Thirty-four  stations  were  counted  in  which  there  was  no 
consistent  time  relationship  between  Hr  and  H.  Only  six 
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examples  were  counted  in  which  and  H  were  strongly  in  phase. 
Figures  20  and  21  show  one  such  station  at  Watheroo,  Australia. 
The  five  remaining  stations  are  Zaimishche-Kasan ,  Russia; 
Kotchlno,  Russia;  Coimbra,  Portugal;  San  Juan,  Puerto  Rico; 
Tananarive,  Madagascar  Island.  Eight  typical  examples  for 
which  there  is  no  consistent  phase  relation  are  Chelyuskin, 
Siberia;  Hongkong,  China;  Vienna,  Austria;  Nantes,  France; 

San  Fernando,  Spain;  Sitka,  Alaska;  Pilar,  Argentina;  Batavia, 
Java.  Eight  typical  examples  for  which  there  is  a  180°  phase 
relation  are  Elisabethville ,  Belgian  Congo;  Yakutsk,  Siberia; 
Toyohara,  Japan;  Tsingtao,  China;  Toungoo ,  Burma;  Uccle (Brussels) 
Belgium;  Meanook,  Canada;  Tuscon,  Arizona. 

The  occurrence  of  this  out-of-phase  relation  is 
sufficiently  striking  to  warrant  further  investigation.  Perhaps 
the  effect  is  pure  within  the  core,  and  is  partially  obscured 
at  the  earth's  surface  by  a  superposition  of  different  sources 
together  with  the  quasi-dispers ive  effect  (Chapter  III). 
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Appendix  A' 

The  following  tables  list  the  algebraic  average  of 
A™  and  for  Vestine*s  X  and  Y  analyses  (cf.  ref.  6,  p.  42), 


hereafter 

m 

referred 

n 

to  as  A™ 

1912.5 

and  B  . 

< 

1922.5 

1932.5 

1942.5 

0 

1 

2470 

2839 

2298 
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0 

2 
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-2815 

-3574 

0 

3 
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1258 
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0 

4 
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0 

5 

-2384 
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0 

6 
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1 

1 

97.5 

367.5 
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1 

2 
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1 

3 
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1 
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4 
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4 

6 
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5 

5 
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5 
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6 

6 
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n 
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1 
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The  following  table  lists  the  time  average  (for  the 
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